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Introduction 


Hole size effect 

Experiments on graphite-epoxy laminated plates containing unloaded small holes show that these 
laminates are notch insensitive. That is, the uniaxial strength of these laminates with small holes 
exceeds the strength predicted by a point stress criterion using the stress concentration factor for 
the in-plane stress field. Laminates containing large holes exhibit notch sensitive behavior and 
consequently their strength is reasonably well predicted by the stress concentration effect. This 
hole size effect is manifested both in tension and in compression. Apparently some mechanism 
must cause in-plane stress relief for laminates containing small holes. 

A mechanism to explain relief of the in-plane stress concentration in the vicinity of the hole for 
laminates in compression was suggested by Mikulas (1980). This mechanism is a geometrically 
nonlinear response that couples the high in-plane stresses with the out-of-plane rotations in the 
vicinity of the hole. A geometrically nonlinear response is associated with reduced growth of the 
in-plane stresses as the applied compressive load is increased since, in effect, the in-plane load is 
directed out-of-plane. Out-of-plane rotations are larger in the vicinity of the hole because of a dif- 
ferential expansion in the thickness of the laminate, or bulging, caused by Poisson’s effect. 

Objective 

The purpose of this research was to study the influence of geometric nonlinearity on the micro- 
mechanical response of a filamentary composite material in the presence of a strain gradient 
caused by a discontinuity such as a hole. A mathematical model was developed at the micro- 
mechanical level to investigate this geometrically nonlinear effect. , 


Mathematical model 

The mathematical model was developed for a solid rectangular parallelepiped, or block, of length 
L, height H, and of a width that is large relative to dimensions L and H. For a cartesian coordi- 
nates ( x,y,z ) with origin at the center of the block, - L/2 <x< L/2 and -H/2 < z < H/2. The block 
is made from continuous fibers parallel to the x-axis that are identical and equally spaced. These 
fibers are encapsulated by a matrix material. The block is subjected to uniform end shortening A/2 
at x = ± L/2, traction-free conditions at z = ± H/2, and body forces are neglected. For a very wide 
block with loads and material properties uniform in the y-direction, plane strain in the x-z plane is 
valid. Let «, v, and w denote displacements in the x-, y-, and z- directions, respectively. Then, 
plane strain implies u = u(x,z), v = 0, and w = w(x,z). Due to symmetry only a quarter of the x-z 
plane is modeled as shown in Fig.l. To simulate an in-plane strain gradient, we prescribe the 
u(x,0) subject to the conditions that u(x,0) is continuous for 0 <x < L/2, u(0,0) = 0, and u(U2,0) = 
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Fig. 1 Solution domain for a fiber reinforced solid block subjected to uniform end 
shortening, traction free conditions on the remaining external surfaces, and sym- 
metry boundary conditions along the coordinate axes. 


-A/2. Within this class of functions, u(x,0) is prescribed such that the magnitude of the normal 

strain e x (x,0) ( = JL [ u (*, 0) ] ) is a local maximum at x = 0, and that this local maximum is 

larger than the average shortening A/L. Prescribing u(x,0) in this manner is the method by which 
an in-plane strain gradient is introduced into the uniform compressive strain field. 


A global and local region are defined for the model as shown in Fig. 1. In the local region the 
response of individual fiber-matrix units is determined by functional degrees of freedom for each 
unit, such that the total number of degrees of freedom depends on the number of fiber-matrix 
units. In the global region the response of individual fiber matrix units are determined by only a 
few functional degrees of freedom defined over the entire global region, and these global degrees 
of freedom are independent of the number of fiber-matrix units. The global region represents a 
smeared model of the individual fiber-matrix units remote from the stress concentration at the ori- 
gin. Details of the fiber matrix units are shown in Fig. 2. 


By approximating the displacements in the thickness coordinate z, the dependence of the field 
equations on the two coordinates x and z is reduced to a set of equations with one-dimensional 
dependence in x. Details of the mathematical development are given in Appendix A. The resulting 
equations are summarized in Eqs. (73) to (77) in Appendix A. There are 8N+8 first order, nonlin- 
ear, ordinary differential equations subject to 8N+8 boundary conditions, where N denotes the 
number of fiber-matrix units in the local region. The numerical solution to this nonlinear two- 
point boundary value problem was attempted by using the computer code called PASVART(Len- 
tini and Pereyra, 1977). Persistent numerical ill-conditioning of the system matrices were encoun- 
tered using this code, so that no acceptable solution could be found. 
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Fig. 2 Configuration and nomenclature of the fiber-matrix units. 
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Appendix A 
Governing Equations 


/* 


Kinematic assumptions 


Local region (n = 1, 2 ,..., N) 

Let tjrand t m denote the thickness of the fiber and thickness of the matrix in a typical fiber-matrix 
unit. The total unit thickness t = tf+ t m , and the volume fractions are Vf = tj/t and V m = t l7 /t. The 
displacements in a typical fiber-matrix unit are assumed to be explicit functions of the thickness 
coordinate z. It is convenient to use dimensionless thickness coordinates in the fiber and matrix. 
For the matrix in the nth fiber-matrix unit, n = 1, 2,..., N, let 


C = (z-z 2 n - 3 )/t„ 

and for the fiber let 


Z2n -3~Z — Z-2n - 2 0^*1 


( 1 ) 


1 „ „ 1 
2 1 2 


^ = (z-z 2n _ l ) /t f z 2n _ 2 — ^ — z 2n _ i ~<ri<^ ( 2 ) 

where z 2n _ l = (z 2n _ 2 + z 2n _ l ) /2. For the nth unit the displacements are assumed to be 

+ C ( 1 - 0 “ 2„-2 O ) 


u (x,z) = " 2 „- 3 -| 0 2 „-3 + ? 


tf 

W 2n-l + 5 0 2n -l- H 2«-3 + 2 0 2 n -3j 


u(x,z) = u 2n _ l -y\tp ln _ l (4) 

w{x,z) = w 2n _ 3 + C[w 2 „_ 1 -w 2n _ 3 ] (5) 

w(x,z) = w 2n _ l ( 6 ) 


There are three functional degrees of freedom for the nth unit. The fiber displacements are U 2 n . 
](x) and W 2 n .j (x), and the matrix displacement is W2«-2W- Fib er rotation is 
dw 2n _ j 

_ j = ^ = w' 2n _ j , where a prime denotes ordinary derivative with respect to x. These 

displacement assumptions given by Eqs. (3-6) satisfy continuity between adjacent units. The 
functional forms for the matrix, Eqs. (3) and (5), are the fewest number of terms in a power series 
in z (or Q which model thickness shearing deformation and thickness stretch. The assumption for 
the fiber displacements in Eqs. (4) and (6) are those of the Euler-Bemoulli beam theory. 

For the first fiber-matrix unit (n =1) symmetry reduces the displacement assumptions in matrix 0 
to 
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( 7 ) 



U(X,Z) = Mj + +C(1 ~C)«0 f-^- 1 

w(x,z) =-w l + 2C ) w l 


Global region (m = N+l, N+2,..., N = M) 

The global region consists of M fiber-matrix units. The thickness of the global region is denoted 
Iiq so that Iiq = Mt. The middle of the global region is at z = Zq = Z2N-1 + h(/2. We define a 
dimensionless thickness coordinate ^ in the global region by 

\ = ( z-z G )/h G ~ (8) 

A separate displacement field is assumed for the global region that is explicit in the thickness 
coordinate. Introducing the five functional degrees of freedom uq(x), \| i x (x), OfxJ, wq(x), and 

\j f z (x), the displacements in the global region are 


u(x,z) =u G + h G tw x + j h G ^ 2 O 
w (x, z) - w G + h G fyy z 

Displacement continuity with the top fiber at each x in the local region requires 

tf /l/j Hr 

U 2N-l~ 2®2N-1 ~ U G~~2^x + 

hr 


w. 


2 AT- 1 


l G 


(9) 


( 10 ) 


(ID 


Point matching the displacements at the top of a typical layer in the global region, using Eqs. (4) 
and (6) with n — > m, to the global field in Eq. (9) gives 


l 2m- 1 ~ 2 0 2m-l “ U G + h G^2m~ l^x + 2^2m- 1® 


( 12 ) 


w 


2m - 1 




where ( = (z-> m _ j - z G ) /h G . Equations (10-12) can be solved for the typical unit fiber 
displacements in terms of the global displacements to get 


t f t f 


l 2m 


- 1 “ U 2 N- 1 2 ^ + ^2/n- 1^ ®2N- 1 + 2 ^ + ^2m- 1^ ®G + 


(13) 
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( 14 ) 


W 2m-1 “ ^ + ^2m-0 W G 2m-\ W 2N-\ 

®2m-l = ^ + ^2m-0^G~^2m-\^2N-\ (15 ^ 


dw r 

where 9 „ = — = w' r . The axial displacement at the center of the matrix in the /nth unit in the 

G dx ° 

global region, obtained from Eq. (3) by letting n -» m and £ = 1/2, is equated to the axial displace- 
ment at the center of the /nth unit from the global field Eq. (9). Doing this we get 


1 


+ ^9- 


-U. 


2 \_ U 2m - 1 + 2 °2 m - 1 + U 2m - 3 2 U 2»j ~ 3 J 


1 


+ 4 M 2m-2 


-2 _ U G + 2 ^2m- 1 + ^2m-3^ ^ + 


T«2„-. + ^-3) 2, » 


(16) 


Equation (16) is not a point matching condition between the /nth fiber-matrix unit and the global 
field since t^/2 * t/2. Using Eqs. (13-15) the final expression for the matrix axial displacement 

u 2m-2 is 


5“2„,-2 = -5 / [(i + 25 2 m .,)e c -25 2 „. 1 9 2K .,]-^® (17) 

Equations (13-15) and (17) relate the fiber and matrix displacements for an individual unit in the 
global region to the global degrees of freedom. 




Strains 


Local region (n = 1, 2,..., N) 

Assuming small strain and moderate rotations in the solution plane, the strain displacement rela- 
tions are 


e 


X 


du 1 ,3w> 2 
dx + 2 ^dx 


du dw 

e z~d~Z y K~Tz + frc 


(18) 


Substituting the displacements in Eqs. (3) and (5) into Eqs. (18) we get the strains in the matrix of 
the nth unit, n > 1, as explicit functions of Theses matrix strains in the nth unit are 


? 2,-2+W-l)e 2 „. 2 +(2C- 1)^.2 

(19) 


(20) 

V, = i- 2 +(2 ? - 1 )i ; 2 »-2 

(21) 
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in which the generalized strain-displacement measures are defined as 


'2n-2 “ 2 


u 2/1- 1 + 2 ® 2n-2 + U 2/1-3 2 ^ 2n ~ 3 + 4 “ 2/i -2 + g [®2n - 1 + ®2/i - 3 ^ ^ 


1 


tr t* 

+ /9’ , + /0' 


2 / 1-2 - 2L 2n-1 2 2n-l 2 " -3 2 2 " _3 j 4 L 2 ” _1 2 n -3 


+ 7[© 2 2 n -,-0 2 2„-^] (23) 


^2n-2 ~ 4“ 2/1-2 + g f®2n- 1 ® 2 /i- 3 ^ 


(24) 


U 2n-\ U 2n-3 , t , Q q v 

Yln-2 ~ ( + 2 1 ^2/i-l +0 2/i-3^ 

m m 


(25) 


U 2n - 2 ^ 1 /Q Q x 

Y2/I-2 “ t + 2^2/i-l ®2/»-3^ 




(26) 


The matrix strains for the first unit (n =1) are determined by substituting the displacements in Eqs. 
(7) into Eqs. (18). These strains have the same distribution in £ as shown in Eqs. (19-21) for the 
nth unit, but the generalized strain-displacement equations are different from those given in Eqs. 
(22-26). For the n =1 unit the generalized strain-displacement equations are 


e 0 = u 'i + { Q \ + \ u 'o 


(27) 


*o = 0 g o = -\ u 'o + I q2 i 


(28) 


2w x 

E = - — 


(29) 


% = 0 T O = -r +0 i 


(30) 


The fiber strains are determined as explicit functions of T| by substituting the displacements in 
Eqs. (4) and (6) into Eqs. (18). For the nth unit the fiber strains are 


v, = ° v , =0 


e = e 

r 2„-t 


(31) 


in which the generalized strain measures for the fiber are defined by 

^2n-l = U 2n - 1 + 2 0 2n - 1 K 2 n - 1 “ ® In - 1 


(32) 


Strains for the global region (m = N+l, N+2,..., N+M) 

The distribution of the strains through the thickness of the mth unit in the global region is assumed 
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to be given by the same expressions for the typical unit in the local region. That is, for the matrix 
in the mth unit of the global region the strain distributions are given by Eqs. (19-21) with n->m. 
For the fiber in the mth unit of the global region the strain distributions are given by Eqs. (31) with 
n -» m. The generalized strain measures in the matrix of the mth unit in the global region are then 
given by Eqs. (22-26) with n — > m, and the generalized strains for the fiber in the mth unit are 
given by Eqs. (32) with n — » m. Substituting Eqs. (13), (14), (15), and (17) into these generalized 
strain measures for the /nth unit, we get the generalized strain measures as functions of the vari- 
ables u' 2 N _ v Q' 2N _ v 9' g , V x > W 2N- i’ and W G- Written in matrix form, the generalized 

strain measures in the mth unit of the global region are 


^2m - 2 

^2m - 2 

^2m - 2 

%m-2 

^2 m - 2 
£ 

Z 2m-2 

^2 m - 1 
K 2m - 1 


0’ 


IQ 


\w 


2 N- 1 
2 N- 1 
0 G 

V* 

O’ 

2 N- 1 




2e c 9 2*-tl 


9 


IN- 1 


(33) 


where the 8 x 7 array of elements C™ and the 8 x 3 array of elements B ™ for the mth unit in the 
global region are listed in Appendix B. 


Equilibrium equations 

The equilibrium equations are obtained from the principle of virtual displacements, or virtual 
work. The external virtual work is zero for the boundary conditions shown in Fig. 1. Thus, equi- 
librium equations are determined from vanishing of the internal virtual work per unit width (y- 
direction); i.e., 8 W int = 0 for every admissible variation of the displacement field, which is given 

by Eqs. (3-6). Let <5 x , <5 z , and x x , denote the stress components for the solution plane. Then 

L 

2 H 

5W /„/ = JJo (° x te X + ° z te z + X xz 5 W dzdx = 0 (34) 

0 

The integral through the thickness is written as a sum of integrals over the individual fiber-matrix 
units. Then the virtual strains, or variation of the strains, are expressed as explicit functions of the 
thickness coordinate using Eqs. (19-21) for the matrix and Eqs. (31) for the fiber. After integrating 
through the thickness of a typical unit the virtual work functional becomes 
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L 

2 N + M 

5W in[ = J ^ {t m [a5e + a5e + a5e + T5Y+x8y+a z 58 z ] 2 ._ 2 + [W5£ + M5 k] 2; _ J dx 
o ' = i 

(35) 

where we defined the matrix stress resultants 


[a, a, a] 2/ _ 2 = J [1, (2£- 1), (2£- l) 2 ] 




= J 


i 5M 


j [1.(2C-D1 v,.5 


*n - Oo ' 


(36) 


and the fiber stress resultants 


1/2 




- 1/2 


(37) 


The variation of the generalized strain-displacement relations, which are obtained from the varia- 
tion of Eqs. (22-30) for the matrix and Eqs. (32) for the fiber, are substituted into Eq. (35). Inte- 
grations by parts with respect to the coordinate x are performed. Vanishing of the functional leads 
to the following equilibrium equations. 


S 2 ’" 2 + ^-l=0 n = 2,3,...,N 


dNo n - 1 

— 2n 1 +x, -f, , = 0 n - 1,2, ...N- 1 

dx 2n 2n ~ 2 


dV 7n . 

-— 2n ~ l + o -a =0 n = 1,2,..., N — 1 

UX *2* - 2 ^2n 


N + M 


5”'' I (c5\.-2 +c 5\.-*)-o 

m = V+ 1 


(38) 

(39) 

(40) 

(41) 


NAG- 1-201 


10 


ly N+M 

aV 2N-\ , ^ . , V r m rr - n 

4* <7 + f / C/rrCJ — U 

Z lN-2 m ~ 6 Z 2m-2 


rfjt 


m = W+l 


( 42 ) 


dS G N + M m _ 

X < C «S»-2 +C ?5%»-2> =0 


tlx 


w = N + 1 


(43) 


dlV, 


N+M 


fa ~ l m X ( C 44 ' C 2 m -2 + C 54 ^ 2 m- 2 ^ “ 0 
m-N+\ 


(44) 


dx 


N+M 


+ t m 


x «v,=» 


m = N + l 

The additional stress resultants appearing in Eqs. (38-45) are defined by 


(45) 


$ 2n-2 ~ 4 ? w( a ^2n-2 « - 1,2, ...,1V 


(46) 


1V 2 „_i — lV 2n _ i + " 2 " [ + O) 2n-2 + 2n^ « — 1,2, ...,A1 1 (47) 

M 2 „-i = M 2 „_, + i(„< / [(o + &) 2 „. 2 + (-5 + 0),.] « = 1.2 W-l «« 


/ '-^ f2 "' 1 +[^-| 8 2 „-,-^- 2 ( ® 2 .-l- 0 2 »- 3 )+i 2 „ ( e 2 „ + ,- 02 .-l ) ] + 


2/j-l 


2^ X 2n-2 + X 2r) + ^2 ( X 2n-2 ^2n^ H- 1,2, ...,1V 1 


N 2 N - 1 “ ^ 2 N- 1 + y ( a + 2 N -2 + 


N + M 


X [C7i" 2 »,-i + ' m (CnS 2 „_ 2 + Cj ] 5 2m . 2 + C”a 2 „. 2 )] 

m = N+ 1 


(49) 


(50) 


1 


M 2 ^_ 1 - M 2N - 1 + + G ) 2W-2 + 


N + M 


X [ C 82^2ffl - 1 + C 12 N 2m - 1 + l m ( C V2 C! 2m - 2 + C 22 (5 2m - 2 + C 32°2m - 2^ 

m = N+l 


(51) 
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2N - 1 


— + ^2N-\^2N-\ ^2N-\^2N-\ ^2N-J + , 2 $G®2N - 1 


^ —N + —S ^ 
7T ; G ,2^0 
V G h G j 


0 G + 0 X 2N-2 + 2'^2//-2 + l m X ^ C 42 X 2m-2 + C 52 X 2m-2^ 

m = N+ 1 

(52) 


N + M 


N + M 


S G “ X ^75^2OT-l + / «(^15 a 2m-2 + ^25®2m-2 +C '35^2m-2) (53) 

m = V+ 1 


V+AT 


W c = X - ! + '», < Cw^m - 2 + - 2 + C 34 & 2,„ - 2> «4) 

m = A/ + 1 


iV+Af 


Mq - X ^83^2m- 1 + ^73^2m- 1 + ^m^l3 a 2m-2 + ^23^2m-2 + ^33^2m-2) ^ 

m = W+l 


dM n 9 2 8 ^ 


V 


G > 


N+ A/ 

( 0 G _0 2iV-I^ + / m X ( C 43"2/n-2" r v '53 k 2m-2 
m = W + 1 

(56) 


m r +r m X ) 

x 2m-2 + c 53 x 2m-2' 


The boundary conditions consistent with virtual work and for the particular problem shown in 
Fig. 1 are 


«2n-2 = 0 

at x = 0, L/2 

n = 

2 , 3, ..., N 

(57) 

(0) = 0 

u 2n -i(L/2) =- 

-A/2 

li 

►— * 
N> 

(58) 

0 2„-l =0 

at x = 0, L/2 

n = 

1 , 2 ,..., A 

(59) 

23 

1 

II 

o 

at x = 0, L/2 

n = 

1 , 2 ,..., A 

(60) 

II 

II 

e 

0 G = ^G - 0 

at x = 

: 0, L/2 

(61) 


At z = 0 the axial displacement u(x,0) is prescribed subject to u(0,0) = 0 and u(L/2, 0) = - A/2. 
Denote the function U(x) as this prescribed displacement function. From the first of Eqs. (7) with 
C, = 1/2, the prescribed displacement relates u\, Gj, and u 0 of the first fiber-matrix unit by 


U(x) = u^x) +t J Q l (x) +i« 0 (x) (62) 

Equation (62) is used to express the matrix displacement function w 0 in terms of the prescribed 
function and the displacement and rotation of the first fiber. For «j and 0j satisfying Eqs. (58) and 
(59), U properly prescribed, matrix displacement uq will vanish at x = 0, L/2. 
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Hooke’s law 


For plane strain of an isotropic material, Hooke’s law for the matrix is 


C , = K m + °z = K m K + Vi 1 = G »J, 


in which 


E (1-0 ) v 

K _ ”» V m } yy _ 

OT ( 1 + x> ) ( 1 - 2d ) m 1 - v 

' m ' v /w y 7, 


2(1+1>J 


and ^ and v m are Young’s modulus and Poisson’s ratio, respectively, for the matrix. 

For the Euler-Bemoulli kinematics of the fiber, Hooke’s law is 

o x = Ej£. x where E f = V f E f + V m E m (65) 

and E { is Young’s modulus for the fiber, V { is the fiber volume fraction, and V m is the matrix vol- 
ume fraction (V f + = 1). The modulus Ej represents the composite modulus of a fiber-matrix 

sheet parallel to the x-y plane of thickness tf. 

Hooke’s law in terms of the stress resultants and generalized strains is obtained as follows. The 
strain distributions given by Eqs. (19-21) for the matrix are substituted into Eqs. (63), and these 
results are in turn substituted into the stress resultant definitions for the matrix, Eqs. (36), to get 


°2i-2 = U - jS/i) [« + 58 + V; 


= - A K m 

4 M 


s o + 2 e o + Vv 


®2(-2 =3^21—2 i>l 

*»->■ 0-I*, 1 ) jr -[5 i+ 3 t+ K e «L. 


T ° 4 G ^° (69) 

^2i-2 = ^njli-2 i > ^ 

hi-2 = ( I -2*n)5 G -Ta-2 '' al < 70 > 

‘- 1 <71> 

The strains given by Eqs. (31) for the fiber are substituted into Hooke’s law given by Eq. (65), and 
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these results in turn are substituted into the resultant definitions in Eqs. (37) to get 

N 2i - 1 = fyfai-l M li-\ = 12^// K 2«-1 ' “ 1 (72 

The task remaining is to substitute the generalized strain-displacement relations for the local 
region (Eqs. (22-30) and (32)) and the global region (Eqs. (33)) into Eqs. (66-72). These results 
are then substituted into the stress resultant definitions in Eqs. (46-56), and into the non-differen- 
tial terms in equilibrium Eqs. (38-45). Although the results of these manipulations are not pre- 
sented here, the form of the stress resultant-displacement relations are presented in the following 
section. 


Summary of mathematical model 


Let the (4N+4) x 1 generalized displacement vector be defined as 

q(x) = [w 0 H 1 9 1 W 1 U 2 U 3 0 3 W 3 M 2iV-2 “2At-l 0 2iV-l w 2^-l ° ^ 0 G W q 




(73) 


and the (4N+1) x 1 generalized force vector as 

_ _ _ ..r 

Q(x) = [s 0 N { Mi V l S 2 N 3 M 3 V 2 ... $2N - 2 ^2 N - 1 1 ^2N-^G N C M G V g\ (74) 

The mathematical model can be posed as a system of 8N + 8 nonlinear, first order, ordinary differ- 
ential equations in the generalized force and displacement vectors subject to boundary conditions 
at x = 0 and x = U2. 

Equilibrium equations: There are 3N+2 equilibrium equations given by Eqs. (38-45). After 
using Hooke’s law and the generalized strain-displacement relations, these equilibrium equations 
are of the form 



(75) 


in which A is a (3N+2) x (4N+4) matrix containing ones and zeros, and is a (3N+2) x 1 vector 
functional of the generalized displacements and their derivatives. 

Generalized force-displacement equations: The generalized force definitions given by the 
4N+4 Eqs. (46-56) can be written in terms of the generalized displacements and their derivatives 
by using Hookes’ law and the generalized strain-displacement relations to eliminate 

ir, cr, a, f, x, a z , N, and M. The result is 4N+4 equations of the type 



(76) 


in which D is a (4N+4) x (4N+4) matrix of constants, and F 2 is a (4N+4) x 1 vector functional of 
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the displacements and their derivatives. The non-zero rows of matrix D result from the definitions 
ofV 2n _i (49), V 2 N-1 (52), and V G (56). 

Kinematic definitions: The rotations are related to the displacements by the equations 
dw i dw r 

— = 9„ , ,and — = 0~. These definitions result in N+l linear equations of the form 

dx 2n ~ l dx G 

J = 7 q (77) 

dx 

in which 7 and 7 are (N+l) x (4N+4) matrices containing ones and zeros. 

Prescribed displacement condition: The axial displacement U(x) is prescribed in Eq.(62) and 
this relates displacements in the first fiber-matrix unit. This single equation completes the set of 
8N+8 equations for as many unknowns. 

Boundary conditions: There 8N+8 boundary conditions given in Eqs. (57-61) for the 8N+8 sys- 
tem of first order differential equations. 


Appendix B 

Matrix elements in Eq. (33) 

The superscript m is dropped for array elements C ™ and B™ in the following formulas for conve- 
nience in writing, but it is implied. 

^11 = 1 ^12 = ^14 = (^ G ^) (1 + + ^ 2 m - 3 ^ 

c 15 =(a£/ 8)(- i + (^ m _,+4 2m . 3 ) 2 ) c 13 = c 16 = c„ = o 

^22 = ^23 = ( l f( 1 + 2£ 2m - 1) ) ^ C 24 = t/2 ^ 

C 25 = (h G / 2) 2 (^2m-l~^L-l) C 21 = C 26 = C 2? = 0 


C 32 = 1 C 33 = (t/2) ( 1 + 2^ 2m _ j) C 35 = t 2 /8 


^31 ~ C 34 ~ ^36 ~ ^37 ~ 0 


C 42 ~ ( t ? /(^ 2 m-l ^2m-3^ + 1 + ^2/n-3^ ^ 

^43 = [ V ( ^>2m - 1 “ ^2m-3^ + { + ^2m- 1 + ^2m-3^ ^ ^ t m 

C 44 = t/t m C 45 = (^c(^ 2 m-l _ ^2m-3^ // (^ r m) 
C4I = ^46 = C47 = 0 


1 


(B-3) 


(B-4) 
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C 52 “ 4 7^^2m-l ^2m-l + ^2m-3 ^53 “ 2 / ( 1 + 2 ^2m- 1^ + ^2m- 1 ^2m- 

m m 

C 53 = < 2 / (20 c„ = C M = C* = C 37 = 0 
(l -t ) 

n _ n — o '^2m-l ^2m-3 _ r — r — r — C -ft 

_C 66 “ C 67 “ 2 7 L 61 “ C 62 “ °63 “ C 64 ~ C 65 u 

l m 

C 7 j = 1 C 72 = -^(1 + 2^ 2m _ j) C 73 = ^(l+2^ 2m _j) 

C 74 = (* c /2) (1 +2^ m .,) c 75 = (*J/8) (4^2m- 1 ~ 1) c 76 = C 77 = < 

C 82 = -2 ^2m-l C 83 = 1+2 ^2 ot-1 C 81 = C 84 = C %5 ~ ^86 = C 87 = 0 

-8ll= (1 + ^>2m- 1 + ^2m-3^ ^12 = + ^2m— 1 + ^2m-3^ ^2m - I + ^2w-3^ 

5 13 = (^2m- 1 + ^2m-3^ 

5 21 = 2 ^2/n- 1 ~ ^2m-3^ ^ + ^2m- 1 + ^2m-3^ 

B 22 ~ ~^2m-l ~^2m-3^ ^ + 2 ^2m - 1 + 2 ^2m - 3^ 5 23 = 2 &2m - 1 “ ^2m - 3^ 

5 3 j = - 8 32 = 833 = (^2m- 1 ~^2m-3^ 2 

8 ;1 = 8j'2 = 8 /3 = 0 / = 4, 5, 6, 8 

871 = (1+2^ 2/jj _ 1 ) 8 7 2 = “2^ 2;n _ 1 (1 +2^ 2to _ 1 ) 8 73 = 4^ 2m _ 1 


3 (B-5) 

(B-6) 

(B-7) 

3 

(B-8) 

(B-9) 

(B- 10) 

CB-11) 

(B-12) 

(B-13) 


NAG-1-201 


16 


